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Pioneering study of Gamow-Teller (GT) and Fermi matrix elements (MEs) using no-core- 
configuration-interaction formalism rooted in multi-reference density functional theory is presented. 

After successful test performed for ®He^®Li /3-decay, the model is applied to compute MEs in the 
sd- and p/-shell T=l/2 mirror nuclei. The calculated GT MEs and the isospin-symmetry-breaking 
corrections to the Eermi branch are found to be in a very good agreement with shell-model predic¬ 
tions in spite of fundamental differences between these models concerning model space, treatment 
of correlations or inclusion of a core. This result indirectly supports the two-body current based 
scenarios behind the quenching of axial-vector coupling constant. 

PACS numbers: 21.10.Hw, 21.60.Jz, 21.30.Fe, 23.40.Hc, 24.80.+y 


The atomic nuclei are unique laboratories to study 
fundamental processes and search for possible signals 
of new physics beyond the Standard Model in ways 
that are complementary or even superior to other sci¬ 
ences. This is due to enhanced sensitivity of specific iso¬ 
topes to fundamental symmetries caused, in particular, 
by intrinsic-symmetry-related or purely accidental (near- 
)degeneracies of nuclear states. For example, parity dou¬ 
blets caused by stable octupole deformation increase sen¬ 
sitivity to search for the violation of CP-symmetry which 
is responsible for matter over anti-matter dominated Uni¬ 
verse [I[. Accidental near-degeneracy of 3/2+ and 5/2+ 
levels in ^^^Th, which are separated only by 7.6±0.5eV, 
opens up a possibility for high-precision measurement 
of the temporal variation of fine-structure constant with 
much higher sensitivity as compared to the atomic tran¬ 
sitions [2(. Last but not least, nuclear physics input is 
critical in an ongoing hunt for a weakly-interacting mas¬ 
sive particle (WIMP), a candidate for dark matter, in 
direct detection experiments measuring the recoil energy 
deposited when WIMP is scattered off the nucleus, see [3| 
and refs, quoted therein. 

Traditionally, the atomic nuclei are used to study the 
weak interaction. A flagship example is the superallowed 
1=0+^/=0+ /3-decay among the members of the iso- 
baric triplets T=l. With small, of order of a percent, 
theoretical corrections accounting for radiative processes 
and isospin symmetry breaking (ISB), these semileptonic 
pure Fermi (vector) decays allow to verify the conserved 
vector current (CVC) hypothesis with a very high pre¬ 
cision. In turn, they provide the most precise values of 
the strength of the weak force, Gp, and of the leading el¬ 
ement, Vud: of the Cabbibo-Kobayashi-Maskawa (CKM) 
matrix, see Q for a recent review. 

The T=l/2 mirror nuclei offer an alternative way to 
test the CVC hypothesis Q. These nuclei decay via the 
mixed Fermi and Gamow-Teller (GT) transitions. Hence, 
apart from the radiative and the ISB theoretical correc¬ 
tions, the final values of Gp and Wd depend on the ratio 


of statistical rate functions for the axial-vector and vec¬ 
tor interactions, /a/Zv, and the ratio of nuclear matrix 
elements p ^ AMgt/VTp where A = QaIqv denotes the 
ratio of axial-vector and vector coupling constants. 

The CVC hypothesis implies that the vector coupling 
constant is a true constant ^v=l- The axial-vector cur¬ 
rent is partially conserved meaning that the coupling con¬ 
stant gets renormalized in nuclear medium. The effective 
axial-vector coupling constant, =(igA^ is quenched 
by an A-dependent factor, q, with respect to its free 
neutron decay value gA ~ —1.2701(25). Quenching fac¬ 
tors deduced from comparisons between the large-scale 
nuclear-shell-mo del (NSM) calculations and experiment 
are: q ^ 0.82, q ^ 0.77 p and q « 0.74 Q in the p-, 
sd-, and p/-shell, respectively. In heavie^ A=100-134, 
nuclei, the average quenching is q ^ 0.48 [Sj. This result 
is consistent, m to a theoretical uncertainty, with the 
result of Ref. [9|. Even stronger quenching, gAA~^-^^, is 
used in the IBM-2 model [l(| . 

The question about physical causes of the quenching 
has no unique answer. The quenching is usually re¬ 
lated to: (i) missing correlations in the wave function, 
(a) truncation of model space or {Hi) to a very fashion¬ 
able nowadays renormalization of the GT operator due to 
the two-body currents Scenarios involving non- 

nucleonic degrees of freedom like NN NA excitations 
are shown to contribute only rather weakly [H. 

Proper understanding of the quenching is essential for 
many branches of modern physics from modeling of as- 
trophysical processes in stars to elusive neutrinoless dou¬ 
ble beta decay which depends on the fourth power of 
and is therefore particularly sensitive to q. In or¬ 
der to address the quenching, it is of paramount impor¬ 
tance to investigate the GT matrix elements (MEs) us¬ 
ing diverse theoretical models. The goal of this work is 
to communicate the pioneering application of the multi¬ 
reference density functional theory (MR-DFT) rooted no- 
core-configuration-interaction (NCCI) approach to study 
/3-decay, with a particular emphasis on the GT process. 
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After a short introduction to the model we shall present 
the numerical results starting with the /3-decay of ^He 
in order to test reliability of the model. Afterwards, the 
model will be applied to the sd- and lower p/-shell T=\j2 
mirror nuclei, where both the GT MEs and the Fermi 
MEs will be computed. 


(non-orthogonal) basis: 

|n; IM- T,) = ^Vi; (1) 

id 

-V V pi I \ 

id K,T>|T^| 


The NCCI models rooted in MR-DFT offer nowadays 
an interesting alternative to the conventional nuclear 
shell model |14l4l6j. Firstly, they are capable of treat¬ 
ing rigorously both the fundamental (spherical, particle 
number) as well as approximate (isospin) symmetries. 
Secondly, by invoking the generator coordinate method 
and/or mixing of discrete (quasi)particle-(quasi)hole (or 
particle-hole) configurations, they allow to incorporate 
important correlations into the nuclear wave function. 
Thirdly, they can be applied to any nucleus irrespec¬ 
tively of A and the neutron and proton number pari¬ 
ties. Moreover, by construction, they are able to capture 
core-polarization effects resulting from a subtle interplay 
between the long-range and short-range nucleon-nucleon 
forces, what is of critical importance for the calculation 
of isospin impurities and ISB corrections 0 - 


The NCCI formalism developed by our group lH, fl^ 
involves the angular momentum and isospin projections 
and subsequent mixing of states having good angular 
momentum and properly treated isospin. It proceeds in 
three distinct steps. First, we compute self-consistently 
a set of k Hartree-Fock (HF) (multi)particle-(multi)hole 
configurations, {\^j)}j=i^ relevant for the problem un¬ 
der study. The Slater determinants {\^j)}j=i are cal¬ 
culated using true Skyrme interaction in order to avoid 
singularities at the next stage, at which we apply the an¬ 
gular momentum and isospin projections to determine, 
for each j, the family of states {\(pj] IMK] TT^)} having 
good isospin TT^, angular momentum IM^ and angular- 
momentum projection on the intrinsic axis K. Subse¬ 
quently, the states {\ipj]IMK]TTz)} are mixed in order 
to account for the K and isospin mixing. This gives a set 
of good angular-momentum states 

for each HF configuration j. The set is non-orthogonal 
and, in general, over complete. In the final step, the states 
are mixed over different configurations by solving the 
Hill-Wheeler-Griffin (HWG) equation, Hu = ENu, with 
the same Skyrme interaction that was used at the HF 
level. The HWG equation is solved in the collective space 
spanned by the natural states corresponding to non-zero 
eigenvalues n of their norm matrix N. The same tech¬ 
nique is used in the code to handle the iC-mixing alone 
as described in detail in Ref. 0 - 


On exit, the NCCI code provides eigenfunctions that 
are labeled by the index n numbering eigenstates in as¬ 
cending order according to their energies and the strictly 
conserved quantum numbers /, M, and Tz = {N — ^)/2. 
The eigenstates can be decomposed in original projected 


where ^ and P^x stand for the isospin and angular- 
momentum projection operators, respectively. This form 
is particularly useful to compute MEs of the GT operator: 


= T(n ; I'M'- /M; T,). (2) 

where ^ expressed by means of 

one-body spherical tensors. The isospin index above is 
fixed /i = ±1 and it determines the overall phase fac¬ 
tor. The matrix element ([ 2 ]) fulfills the Wigner-Eckart 
theorem: 

= -;jA= VVi. (n', J'llO^lln, J) (3) 


where stands for the Clebsch-Gordan coefficient. 

The reduced matrix element equals: 




E Jn'd'M'-p)* An-dM-T^) 
Ji'j'K'T JijKT 


ijKT i'j'K'T' 


V 2 /' +1 a 


T'T' 


-77,177 


rii'K' 






The integral, J, runs over the beta Euler angle in 
isospace, /3t, and the Euler angles in space ft = (a, /3, 7 ): 


"^vPdd' 


2 T 


1 r 
Jo 


dpT sin (3t 


2/-M 

Stt^ 


/ 




(5) 


where d^, and ^ Wigner functions and cpj 

denotes Slater determinant rotated in space and isospace. 
Mean-field (MF) matrix elements in Eq. (|5]) can be ex¬ 
pressed by means of transition densities. The formu¬ 
las are somewhat lengthy and will be published in our 
forthcoming paper. All integrals appearing above can be 
calculated exactly by applying appropriate quadratures, 
see 0 for further details. 

The Wigner-Eckart relation m implies that the total 
probability of decay summed up over the components u 
of the operator and over polarizations of the final 
state M' is: 




din'TWOdnJd _ 2 \Mgt\ 

— 9 A or , 1 = 9 A ~ 


2 / + 1 


2 / + 1 


(6) 


The calculations discussed below were done using a 
new unpublished version of the HFODD solver [l 8 |, |l9| , 
which was equipped with the NCCI module. In order to 
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track mean-field (MF) configurations and facilitate con¬ 
vergence properties at the MF level all reference states 
used in the NCCI calculations has been self-consistently 
calculated assuming parity and signature symmetries. In 
the calculations we have used a basis consisting either 
N=10 or 12 spherical harmonic oscillator (HO) shells. 
The calculations has been performed using either the SV 
Skyrme force of Ref. [ 2 ^ or a variant of this force, dubbed 
SVso, having 20% stronger spin-orbit interaction. The 
latter force was introduced to improve slightly the single¬ 
particle (s.p.) properties of the SV, in particular, by 
shifting the d ^/2 subshell with respect to the Si /2 sub¬ 
shell. The near-degeneracy of these subshells in the SV 
energy density functional (EDF) causes strong mixing 
and, in turn, leads to unphysically large ISB corrections 
in the superallowed 0+—>0+ /3-decay in H=38 isospin 
triplet [l7|. The SVso force does not cure the problem 
of incorrectly placed s.p. levels but makes the spectrum 
slightly more realistic. Indeed, in ^^Ca, the 51/2 and d ^/2 
neutron subshells are separated by 0.22MeV in the SV 
EDF and l.lSMeV in the SVso EDF, respectively. For 
comparison, the experimental splitting is 2.55MeV 0. 



FIG. 1. (Color online) GT ME for beta decay of ®He in func¬ 
tion of a number of configurations taken in the NCCI calcula¬ 
tions in ®Li and ®He. Solid line shows the experimental value 
of Ref. [ 13 . 

Recently, Knecht et al. [HI performed high-precision 
measurement of the 0+^1+ beta decay of ^He, which 
proceeds exclusively to the ground state (GS) of ^Li 
and determined the corresponding GT matrix element 
|Mgt 1=2.1645(43) assuming ^a=-1-2701(25). This is an 
excellent test case for our model mainly because of a lim¬ 
ited number of ph configurations that can contribute in 
these p-shell nuclei. The results of the NCCI calculations 
performed for this transition are depicted in Eig. [H The 
figure shows the calculated GT ME versus a number of 
configurations taken in the mixing. The very left point 
corresponds to a situation, where no mixing was per¬ 
formed in neither of the nuclei. In this limit, called here¬ 
after the MR-DET limit, the HE reference states were 
selected based ultimately on the energy criterion. Note, 
that already in this limit the calculated ME is in fair 
agreement with the empirical value underestimating it 
by ^7%. Next, keeping the wave function of ^He fixed, 
we have attempted to correlate the wave function of ^Li 


by admixing 1+ states projected from the lowest ph con¬ 
figuration (second point) and from the first two lowest 
ph configurations (third point). This caused an increase 
of the ME to 2.208 and 2.223, respectively, i.e. circa 3% 
above the experiment, see Eig.[TJ At this point we freeze 
the wave function of ^Li and attempt to correlate the 
wave function of ^He (last two points). This weakly in¬ 
fluences the ME giving eventually 2.238. The test shows 
that the model is capable to capture main features of the 
wave functions that are important for reliable reproduc¬ 
tion of the GT ME and provides stable predictions in 
function of a number of admixed configurations. 



FIG. 2. (Golor online) Theoretical binding energies of T=l/2 
mirror nuclei calculated using the NGGI framework. The re¬ 
sults are shown relative to the experimental data. 

Encouraged by the result obtained for the ^He decay, 
we have performed systematic study of the GT and Eermi 
MEs for the GS^GS transitions in the T=l/2 mirror 
nuclei covering the sd- and lower p/-shell nuclei from 
A = 17 till 55. All results shown below were obtained 
using the SVso EDE. This functional, apart from slightly 
more realistic s.p. levels, is also superior in reproducing 
binding energies (BE) in comparison to the SV EDE. The 
ability to reproduce masses is considered to be one of the 
most important signatures of the quality of DET-rooted 
models. The calculated BE relative to empirical results 
are depicted in Eig. O Although the theory tends to 
overbind the lightest species and underbind the heavier, 
the overall agreement is at a quite impressive level of 
±1%. It is better almost by a factor of two than the 
level of agreement obtained for the SV interaction. 

It appears also that the SVso has reasonable spectro¬ 
scopic properties. The theory is able to reproduce the GS 
spins already at the MR-DET level with the exception of 
A = 19 case, where the model predicts I = 5/2+ instead 
of /= 1/2+ to be the GS spin. The energy spectra are, 
in general, in fair agreement with data. Eor the sake 
of illustration Eig. [3] shows theoretical and experimen¬ 
tal 1=3/2+ and 5/2+ states in the lower sd-shell nuclei. 
Similar agreement is obtained for heavier nuclei. 

Eig. [4] shows the GT MEs, I^AAfcxI, calculated using 
the MR-DET and the NCCI models. The NCCI calcula¬ 
tions involve typically four-five low-lying ph ME config¬ 
urations. The results of both approaches are strikingly 
similar except for the ^^V^^^Ti transition. Both mod- 
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FIG. 3. (Color online) Excitation energies of the lowest 3/2^ 
and 5/2^ states in sd-shell T=l/2, Tz—^j^ nuclei ranging 
from A=19 till 33. The calculated (experimental) levels are 
marked by thick solid (dashed) lines, respectively. The calcu¬ 
lated states come from MR-DFT. Theoretical energies have 
been normalized to the experimental ground states. 


els systematically overestimate the data beside the nu¬ 
clei ranging from ^4=29 to 35. Strong suppression of the 
GT MEs in this mass range is related to aforementioned 
clustering of the Si /2 and d ^/2 subshells in ME calcu¬ 
lations. Proximity of these two subshells causes strong 
mixing, which has destructive impact on the calculated 
MEs. Comparison of the MR-DFT results obtained us¬ 
ing the SV and SVso EDFs supports this conclusion. In¬ 
deed, the MEs calculated using these two functionals are 
almost identical everywhere except for the mass region 
discussed above. 


most perfectly with the NSM quenching in spite of nu¬ 
merous differences between the models. In particular, our 
model: (z) includes the core and takes into account core 
polarization effects, (ii) accounts for correlations in a dif¬ 
ferent, more schematic way, than the shell-model, (Hi) 
uses functionals, which were not optimized for the NCCI 
calculations or {iv) uses completely different model space. 
In order to address the latter point, we have performed 
additional set of MR-DET calculations using larger basis, 
consisting A^=12 spherical HO shells. We found that the 
increase of the basis size has almost no impact on the 
calculated MEs. Note also, a systematic difference be¬ 
tween the shell-model and NCCI results in the heaviest 
calculated nuclei. The origin of this difference requires 
deeper study. 
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FIG. 4. (Golor online) Gamow-Teller matrix elements calcu¬ 
lated using the MR-DFT (triangles) and NGGI (diamonds) 
approaches in comparison with experimental data (dots) 
taken from Ref. 0 (sd-shell) and Refs. 0, [IS] (p/-shell). 

Fig. [5] shows the NCCI results in comparison to the 
NSM calculations of Ref. @ (sd-shell) and Refs. 0, [HI 
(p/-shell). The two sets of calculations are very consis¬ 
tent with each other. Indeed, linear fit to the MR-DFT 
(NCCI) GT MEs gives g=0.77 (0.78) in the sd-shell (ex¬ 
cluding problematic A = 31 — 35 cases) and ^^=0.75 (0.69) 
in the lower p/-shell, respectively. These values agree al- 


FIG. 5. GT MEs calculated using the NGGI model (trian¬ 
gles) in comparison to the NSM results of Ref. 0] (sd-shell) 
and Refs. 0 , [2^ (p/-shell). 

The Ikeda sum rule is an important indicator of the 
quality of theoretical models. Eor the T=l/2 mirrors 
it takes particularly simple form: H(0+; n, / ^ 

n', /') = 3. Systematic study of the Ikeda sum rule with 
the present formalism involving the isospin and angular 
momentum projections is CPU expensive. Hence, it was 
limited here to one of the simplest cases of A=39 nuclei. 
In this case, inclusion of all possible ph excitations within 
the 5(i-shell exhausts 99% of the sum rule as illustrated 
in Eig. [6l It is worth mentioning that inclusion of ph 
excitations between the spin-orbit partners 6 ^ 5 / 2 ^ 6^372 is 
crucial for the sum rule. More systematic study of the 
sum rules will be done with the variant involving only 
the angular momentum projection. 

The use of NCCI approach involving both the isospin 
and angular momentum projected states is absolutely 
necessary to study Eermi transitions and, in particular, 
to extract the ISB corrections to the Eermi branch of /3- 
decay. Let us recall that the ISB corrections are needed 
to study the CVC hypothesis and the CKM matrix via 
the transitions in the mirrors, see Ref. The results 
obtained in this study are collected in Tab. [H It is bene- 
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FIG. 6. (Color online) Ikeda sum rule for the GS of ^^Ca. 
Thick horizontal lines indicate the theoretical (solid) and ex¬ 
perimental (dashed) total BE of the GS in ^^Ca and the lowest 
1/2^, 3/2^, and 5/2^ states in The numbers over the 
arrows indicate the total calculated GT strength for each I. 
Multiple arrows and shadowing indicate that the strength is 
distributed over several states. 


TABLE I. Theoretical ISB corrections, (in %) 

adopted in this work. Eor the sake of comparison the table 
contains also the NSM results, of Ref. 


A 

r(NCCl) 

"c 

.(NSM) 

'^C 

A 

r(NCCl) 

'^C 

.(NSM) 

"C 

17 

0.166(17) 

0.585(27) 

37 

0.907(91) 

0.734(61) 

19 

0.339(34) 

0.415(39) 

39 

0.318(32) 

0.855(81) 

21 

0.300(30) 

0.348(27) 

41 

0.426(43) 

0.821(63) 

23 

0.316(32) 

0.293(22) 

43 

0.690(69) 

0.50(10) 

25 

0.413(41) 

0.461(47) 

45 

0.589(59) 

0.87(12) 

27 

0.439(44) 

0.312(34) 

47 

0.673(67) 


29 

0.520(52) 

0.976(53) 

49 

0.646(65) 


31 

0.585(59) 

0.715(36) 

51 

0.714(71) 


33 

0.705(71) 

0.865(59) 

53 

0.898(90) 


35 

0.366(37) 

0.493(46) 

55 

0.620(62) 



ficial to see that our corrections are very consistent with 
the NSM results of Ref. [^ . 

In summary, we have presented a systematic study of 
the GS^GS GT and Fermi MEs in T=l/2 mirror nuclei 
using, for the first time, the NCCI approach based on the 
isospin and angular momentum projected MR-DFT for¬ 
malism. The framework is universal and can be applied 
to any nucleus irrespectively on its mass and the proton 
and neutron number parities. It can be also improved 
and optimized in many different ways, in particular, con¬ 
cerning the tensor force which is known to have an impact 
on the shell structure mil and y5-decay [^ . 

In the present implementation with the SV or SVso 
EDFs the calculated GT MEs systematically overesti¬ 
mate experimental data for the free neutron strength of 
the axial current. The level of disagreement is found to 
be very similar to the one obtained using large scale shell 
model in spite of the fundamental differences between the 
two approaches in handling the core and the core polar¬ 


ization effects, the correlations or the basis truncation. 
It strongly suggests, that the mechanism of in-medium 
renormalization of the axial strength may indeed be re¬ 
lated to the two-body currents. This conjecture requires 
further studies. 

Eventually, we have also calculated the ISB corrections 
to the Fermi decay branch in the T=l/2 mirror nuclei. 
The corrections turn out to be in a very good agreement 
with the NSM calculations. 
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